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ABSTRACT 

We combined two-dimensional kinematic and morphology information on the Ha emis- 
sion, obtained using near-infrared integral field spectroscopy, with broad-band pho- 
tometry to investigate the dynamical structure and the physical properties of a sample 
of ten late-type galaxies at 1.0<z<1.5. Their star formation rate ranges from ~4 to 
~400 Mq yr~^with a mean value of ~80 yr~^. We found that three of these objects 
are undergoing a strong burst of star formation. The sample displays a range of kine- 
matical types which include one merger, one face-on galaxy, and eight objects showing 
evidence of rotation. Among these eight objects, half are rotation-dominated galax- 
ies, while the rest are dispersion-dominated. We found also that two galaxies out of 
the rotation-dominated galaxies are pure rotationally supported disks. They achieve a 
maximum velocity of ~ 180-290 km s~^ within ~ 0.5-1 kpc, similar to local spirals with 
thin disks. Regarding the perturbed rotation and the dispersion-dominated galaxies, 
they display a plateau velocity range of 105-257 km s~^, which is certainly underesti- 
mated due to beam smearing. However, their plateau radii (4.5-10.8 kpc) derived from 
our rotating disk model are significantly higher than those derived for pure rotating 
disks and local spiral galaxies. The galaxies of our sample have relatively young stellar 
populations (< 1.5 Gyr) and possess a range of stellar mass of 0.6-5 xlO^^M©. In 
addition, most of them have not yet converted the majority of their gas into stars (six 
galaxies have their gas fraction >50 per cent). Therefore, those of them which already 
have a stable disk will probably have their final stellar mass similar to the present-day 
spirals, to which these rotating systems can be seen as precursors. We conclude our 
study by investigating the stellar mass Tully-Fisher relation at 1.2<z<1.5. 

Key words: galaxies: evolution - galaxies: formation - galaxies: kinematics and 
dynamics - galaxies: spiral - galaxies: stellar content - galaxies: starburst 



1 INTRODUCTION 

In the last four years, observations by using integral field 
spectroscopy (IFS) of individual star-forming galaxies have 
revealed that galaxies at z>1.5 show a large variety of 
kinematic and dynamical properties. Recent studies at z ~ 
1.5 - 1.6 dWright et aLl[2QQ7l [20091 [Bournaud eTaLl[2QQ8 > 
have confirmed that in addition of the evidence of organized 
rotation non-negligible random motions are also detected at 
this redshift. In addition, it appears that large and mas- 
sive disks with strong star formation already exist at z > 1.5 



Based on data obtained with the European Southern Observa- 
tory Very Large Telescope, Paranal, Chile, programs 075.A-0318 
and 078.A-0177. 
t mrlb@astro.ox.ac.uk 



dForster Schreiber et al.||2006| IGenzel et al.||2006| 12008). 



However more numerous studies at z ~ 2 — 3 (Genzel et al 



2QQ6l|Forster Schreiber et al.'2006' "Law et al."2007"Bouche 
et al.||2007| IGenzel et al. 2008, Nesvadba et al. 2008, van 



Starkenburg et al.|[2008| |Law et al.||2009| ) have found an 
increase in non-circular motions relatives to lower redshift 
samples. Indeed most of these high-redshift galaxies have 
high velocity dispersions. Even when a large-scale velocity 
gradient throughout the galaxy is detected, the value of the 
ratio of V/a is relatively low (< 10) suggesting that it is 
unlikely for these galaxies to have a dynamically cold ro- 
tating disk of ionised gas similar to the local spiral galaxies 
{V/a - 15 - 20 fPib et al.||2QQ6|)). Most of these 'heated 



disks' appear also to be extremely rich in gas with evi- 
dence of high turbulent star formation(Lemoine-Busserolle 
et al. 2009). Another significant difference in the proper- 
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Table 1. Table of observations 



Galaxy 


VVDS-ID 


RA (J2000) 


DEC (J2000) 


z(l) 


Field (2) 


Grating 


Exp. Time 


Seeing (3) 


Run ID 


VVDS-1235 


020461235 


02:26:47.110 


-04:23:55.71 


1.0351 


VVDS-O2I1 




J 


l.llh 


o'.'m 


078.A-0177(A) 


VVDS-2331 


020182331 


02:26:44.260 


-04:35:51.89 


1.2286 


VVDS-O2I1 




H 


31i 


o'.'m 


078.A-0177(A) 


VVDS-6913 


220596913 


22:14:29.184 


+00:22:18.89 


1.2667 


VVDS-22I1 




H 


I.75I1 


O'.'Al 


075.A-0318(A) 


VVDS-5726 


220015726 


22:15:42.455 


+00:29:03.59 


1.3091 


VVDS-22I1 




H 


2I1 


0^.55 


075.A-0318(A) 


VVDS-4252 


220014252 


22:17:45.690 


+00:28:39.47 


1.3097 


VVDS-22I1 




H 


2I1 


0^.61 


075.A-0318(A) 


VVDS-4103 


220544103 


22:15:25.708 


+00:06:39.53 


1.3970 


VVDS-22I1 


J 


& H 


Ih (J); Ih (H) 


0'/.69 


075.A-0318(A) 


VVDS-4167 


220584167 


22:15:23.038 


+00:18:47.01 


1.4637 


VVDS-22I1 


J 


& H 


Ih (J); I.75I1 (H) 


0.^77 


075.A-0318(A) 


VVDS-7106 


020147106 


02:26:45.386 


-04:40:47.39 


1.5174 


VVDS-O2I1 




H 


2I1 


0^.89 


075.A-0318(A) 


VVDS-6027 


020116027 


02:25:51.133 


-04:45:04.48 


1.5259 


VVDS-O2I1 




H 


1.6711 


0'/.65 


075.A-0318(A) 


VVDS-1328 


020261328 


02:27:11.049 


-04:25:31.60 


1.5291 


VVDS-O2I1 




H 


Ih 


0^.61 


075.A-0318(A) 



The columns are as follows: (1) redshift estimated from optical spectrum obtained with VIMOS, (2) VVDS-22h wide field 

(17.5 < Iab < 22.5) and VVDS-02h deep field (17.5 < Iab < 24.0), (3) median seeing estimated from the PSF stars (one taken per 

hour of observation). 



ties of some of these high-redshift objects in comparison 
to the local Universe is the increase in the irregular and 
asymmetric shape of these galaxies. Although some of these 
irregular objects may be associated with mergers, the mor- 
phology and kinematics of the majority of them are incom- 



patible with being ongoing mergers (Shapiro et al. 2008) 



and suggest that they are rotationally s upported ([Genzel 



e^al. 2006 ' Forster Schreiber et al.||2QQ6l [Elmegreen et al 



2007 Genzel et al. 



2008). These so-called 'chain-galaxies' 



have the appearance of giant highly clumpy dis ks. The ob- 
jects can be kpc wide and as massive as 10^ Mq (Elmegreen 



k Elmegreen|2005 Elmegreen et al.|20Q9 ). Elmegreen et al 



(2005) was the first to suggest that they could be the pro- 
genitor of z ~ 1 spirals. Although the origin of the 'heated 
disks' remains uncertain and highly debated, recent theoret- 
ical studies ( Dekel et al.||2009a|b ) propose a bimodality in 
galaxy type by ^ ~ 3 with clumpy star-forming disks and 
spheroid-dominated galaxies with low star formation rates 
(SFRs). At z < 1, the disks should be stabilized by the dom- 
inant stellar disks and bulges, suggested that the redshift 
region of 1 < z < 2 is a crucial step in the formation of 
massive disk galaxies. Nowadays theoretical studies try to 
understand the formation and evolution of extended turbu- 
lent rotating disks, and observational studies, as the ones 
presented in this paper, probing the kinematics properties 
(like observed morphology, dynamical mass, etc), the stellar 
population properties (age of the population, stellar masses, 
etc), star forming rates, and the properties of the ionized gas 
(mass de gas and gas fraction) of the 1 < 2; < 2 disks are 
mandatory for setting constrains on galaxy formation and 
evolution. 

At lower redshifts, the velocity fields of spiral/rotating 
disk galaxies have been used to place important constraints 



on total masses and hence on dark matter halo masses ( Con- 



selice et al. 2005). These rotating disks produce a Tully- 
Fisher relationship (the scaling law between total dynamical 
mass and the luminous stellar mass) which has apparently 
not evolved in slope and scatter since z ^ 0.6 (Puech et al. 



2006 2008). Therefore the large scatter found in previously 



reported Tully-Fisher relationships at moderate redshifts is 
produced by galaxies with perturbed rotation or complex 
kinematics. The importance of the role of the non-ordered 
motion through the gas velocity dispersion for galaxies show- 
ing a perturbed rotation was investigated by |Weiner et al.| 
(2006 ) and Kassin et al. (2007). They defined a new tracer of 



galaxy-dark halo potential, which combines dynamical sup- 
port from the rotation with that of the non-ordered mo- 
tion. The new stellar mass Tully-Fisher relation shows no 
detectable evolution neither in slope or of its intercept up to 
Z'^1.2 ( Kassin et al.||200"7 ). Two-dimensional velocity fields 
allow rotation curves to be deduced in a more robust man- 
ner than slit spectra ( Weiner et al.|[2006 ). IFS allows pure 
rotationally-supported disks to be distinguished from other 
dynamically-disturbed galaxies, which include minor or ma- 
jor mergers, merger remnants and/or inflow/outflows. 



In this paper, we explore the kinematic properties of 
the ionized gas and star formation rates of a sample of ten 
galaxies at 1 < 2; < 1.5 selected in the VIMOS VLT Deep 
Survey (VVDS), using integral field unit (IFU) H-hand and 
J-band spectroscopy with VLT/SINFONI. The results pre- 
sented here are part of a study to investigate kinematics 
and physical properties of unbiased samples representative 
of the global intermediate and high- 2; population. This is the 
second paper of three articles based on data obtained dur- 
ing two observing runs with SINFONI. From the two other 
companion papers, one ( Queyrel et aLllToO Q*) diskusses the 
chemical properties of the intermediate redshift sample pre- 
sented here and the other presents the results found on the 
z ^ 3 galaxy sample ( Lemoine-Busserolle et al.||2009 ). 



In Section[2]we describe our sample, observational strat- 
egy and the data reduction techniques. In Section [S] we ad- 
dress the properties of the nebular emission and also the 
nature of the stellar population of the galaxies, using broad- 
band photometry to constrain properties such as stellar 
mass, age and SFRs. In Section [4] we investigate the kine- 
matic structure and dynamical properties inferred from IFS. 
In Section [S] we explore the nature of these star-forming 
galaxies and in Section |6] we diskuss various relations be- 
tween the properties of the gas and those of the stellar pop- 
ulation and compare our findings with previous studies of 
the kinematics of intermediate and high-redshift galaxies. 
Finally in Section [T] we summarize our results and diskuss 
their implications. 



We assume a cosmology with Qo = 0.3,A = 0.7 and Ho 
= 70 km s~^ Mpc~^throughout, and all magnitudes are on 
the AB system ( |Oke k Guiiii||1983[ ). 
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Table 2. Nebular emission properties 



Galaxy 




F(Ha) (1) 


L(Ha) (2) 


(3) 


E(B-V)^as (4) 


SFRcorr (5) 


c^iD (6) 


Mgas (7) 


m(8) 


VVDS-1235 


1.0352 


20.1 ±2 


11.4 ± 1 


5.3 ±0.5 


0.17 ±0.02 


8.8 ± 1 


74 ±9 


8.4 ±0.7 


0.58 ±0.08 


VVDS-2331 


1.2286 


45 ±8 


39.1 ±7 


18 ±3 


0.28 ±0.03 


42 ±8 


94 ±21 


31 ±4 


0.74 ±0.06 


VVDS-6913 


1.2660 


29.7 ±3 


27.8 ±3 


12.8 ± 1 


0.09 ±0.02 


17 ±2 


113 ± 18 


22 ±2 


0.34 ±0.09 


VVDS-5726 


1.2927 


23.4 ±0.3 


23.8 ±0.3 


10.9 ±0.1 


0.24 ±0.02 


23 ±2 


128 ±2 


19 ± 1 


0.28 ±0.07 


VVDS-4252 


1.3099 


201 ± 54 


205 ± 55 


94 ±25 


0.36 ±0.04 


280 ± 83 


129 ± 40 


143 ± 30 


0.87 ±0.04 


VVDS-4103 


1.3966 


208 ± 49 


248 ± 59 


114 ±27 


0.42 ± 0.04 


412 ± 108 


107 ±57 


208 ± 39 


0.93 ±0.02 


VVDS-4167 


1.4656 


33.3 ± 1 


44.6 ± 1 


20.5 ±0.6 


0.10 ±0.01 


28 ± 1 


141 ±6 


34.1 ±0.8 


0.58 ±0.08 


VVDS-7106 


1.5194 


21.2 ±2 


31.1 ±3 


14.3 ± 1 


0.05 ±0.01 


16 ±2 


82 ±7 


16 ± 1 


0.67 ±0.06 


VVDS-6027 


1.5301 


5.63 ±0.6 


8.4 ±0.9 


3.85 ±0.4 


0.00 ±0.02 


3.9 ±0.4 


42 ±7 


5.5 ±0.4 


0.33 ±0.08 


VVDS-1328 


1.5289 


5.48 ±0.6 


8.2 ±0.9 


3.76 ±0.4 


0.01 ±0.02 


3.9 ±0.5 


78± 11 


5.7 ±0.5 


0.34 ±0.12 


le columns are 


as follows: (1) Flux(xl0 ergs s 


^ cm (2) Luminosity ( X 10+^-*^ ergs s 


(3) raw SFR derived from Ha-line 



(M0 yr (4) reddening suffered by the ionized gas, (5) dereddened SFR, (6) global velocity dispersion derived from Ha-line (km/s), 

(7) total mass of ionized gas (10^ Af©), (8) gas mass fraction. 



2 DATA AND OBSERVATIONS 



3 SPECTROPHOTOMETRY PROPERTIES 



2.1 Sample Selection 

The 10 galaxies, presented in this paper, were selected in 
the VV DS-Q2h "deep" field (17.5 < Iab < 24.0) of the 
VVDS dLe Fevre et"aLl[2QQ5]) and in the VVDS- 22h "wide" 
field (17.5 < Iab < 22.5; |Garilli et aLl[2QQ8t. The ob- 



jects of this sample have spectroscopic redshifts, derived 
from VLT/VIMOS spectroscopy, between z ^ 1.03 and 
z ^ 1.53 (see Table[T]). Their redshifts are such that the 
expected wavelengths of the rest-optical Ha and [N II] emis- 
sion lines would be clear of bright OH sky lines in the near- 
infrared (NIR). These targets were mainly selected on the 
basis of their measured intensity of [O ii] AA3726,3728 emis- 
sion lines in the VIMOS spectra. We chose objects showing 
the strongest [OII]3727 emission lines (rest-frame equivalent 
width (EW) > 50 A and flux > 5 x 10"^'^ ergs"^cm-2) and 
having a photometric spectral energy distribution (SED) 
corresponding to late- type star-forming galaxies. 



2.2 SINFONI observations 

The NIR spectroscopic observations were acquired with the 
3D spectrograph SINFONI at ESO-VLT during two 4-nights 
runs, on 2005 September 5-8 (ESO run 75.A-0318) and 2006 
November 12-15 (ESO run 78.A-0177). SINFONI was used 
in its seeing-limited mode, with the 0.25^^ pixel scale leading 
to a field-of-view of 8x8 arcsec^^ and the H grism (1.447 
- 1.847 /xm) providing a spectral resolution R ^ 4000. One 
galaxy, VVDS-1235, has been observed with the J grism 
(1.098 - 1.399 /xm), where the detection of the Ha emission 
line is expected. Three galaxies (VVDS-1235, VVDS-4103 
and VVDS-4167) have also been observed with the J grism 
(1.098 - 1.399 /im). Conditions were photometric and the 
median seeing for each objects, estimated on 2D imaging of 
point spread function (PSF) stars taken each hour of object 
observation, is indicated in Table[l] Observing strategy and 
data reduction have been described in ILemoine-BusseroTIel 



3.1 Nebular integrated emission properties 

3.1.1 One- dimensional spectral properties of Ha emission 

Integrated rest-frame optical one-dimensional spectra were 
extracted from the SINFONI data cube for each ob- 
ject. From these ID spectra, we measured line fluxes and 
linewidths of the integrated Ha line emission (see Figjl]). 
Using the IRAF task 'splot', we derived the systemic red- 
shift ZHa, the total Ha flux and the global velocity dispersion 
aiD (corrected for the instrumental resolution of ^6.5 A in 
the J band and --6.8 A in the H band). Table [2] list the 
properties of the nebular emission of the 10 galaxies of our 
sample. It is interesting to note that all our targets, except 
VVDS-6027, have a value of aiD > 70 km s"^ 



3.1.2 Star Formation rates from Ha emission 

The Ha nebular recombination line is a direct probe of the 
young, massive stellar population and therefore provides 
a nearly instantaneous (i.e. averaged on the last ten mil- 
lion years) measure of the SFR. The Ha line is moreover 
not strongly affected by dust extinction. We have calcu- 
lated SFRhcx following the 'Kennicutt (1998) calibration, 
re-normalized to Chabrier (2003) IMF: 



SFR^o.{Mq yr"^) = 4.6 x 10"^^ L{Ha) (erg s"^). 



(1) 



Table [2] list the star formation rates deduced from the 
Ha emission for the 10 galaxies of our sample. We found 
that the non-dust-corrected SFRua is between ^4 — 114 
M02/r~^, with seven of our galaxies with a SFRua > 10 
Mo^/r"^. These high SFRs are expected, taking into account 



our selection criteria for the [OH] flux (see section 2.1 ) 



3.2 Stellar population properties 

3.2.1 Spectral Energy Distribution and best- fit models 

We derived the spectral energy distributions (SED) of our 
10 galaxies from broad-band photometry. The magnitudes in 
BVRI filters obtained with the CFH12k ca mera on CFHT 
dMcCracken et al.|2003||Le Fevre et al.|2004 ) were extracted 



et al. (2009) 



from the VVDS catalog http : / / cencosw . oamp . f r| ( lovino 
et al.|20()5 Temporin et al.|2008 ), and magnitudes in ugriz 
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VVDS-020461235 



:1.0351 



= 1.0352 



VVDS-020182331 



z^ojj]- 1.2286 Zh«= 1.2286 




3700 3800 6550 
wavelength (A) 



VVDS-220596913 



6600 
a. 2660 



3700 3800 6600 
wavelength (A) 



VVDS-220015726 



6700 
= 1.2927 




- 



3700 3800 6600 
wavelength (A) 



VVDS-220014252 



a. 3097 



6700 
= 1.3099 



3700 3800 6600 
wavelength (A) 
VVDS-220544103 Zr.,n=1.3970 z„ 



6700 
= 1.3966 




3700 3800 6600 
wavelength (A) 
VVDS-220584167 z, 



6700 

= 1.4656 



3700 3800 6600 
wavelength (A) 
VVDS-020147106 Zr.,n=1.5174 z„ 



6700 

= 1.5194 




6600 



3700 3800 

wavelength (A) 
VVDS-0201 16027 Zr.„-,= 1 .5259 



6700 
= 1.5301 



6600 



3700 3800 

wavelength (A) 
VVDS-02026 1 328 Zr.„n= 1.5291 



6700 
= 1.5289 




3700 3800 6600 
wavelength (A) 



6700 



3700 3800 6600 
wavelength (A) 



6700 



Figure 1. One-dimensional rest-frame spectra for the 10 galaxies of our sample obtained with VIMOS (left) and with SINFONI (right), 
in optical and in i^'-band (or J-band for VVDS-1235) respectively. The position of the fohowing hnes are shown by a vertical line, 
from left to right: [OII]A3727, Ha, [NII]A6584, [SII]A6717, and [SII]A6731. The redshifts inferred from [OII]A3727 and from Ha hnes are 

inrliVatprI in tViP titlp nf parVi r>lntQ 
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Table 3. Photometric data and Spectral energy distributions 



Galaxy Bab ^ ab Rab Iab ^ab gAB i'ab Ub ^ab Jab ^ab 



VVDS-1235 


24.1 


± 


0.06 


23.7 


± 


0.05 


23.4 


± 


0.05 


22.6 


± 


0.05 


24.1 


± 


0.03 


23.8 


± 


0.02 


23.5 


± 


0.02 


22.8 


± 


0.01 


22.7 


± 


0.04 


22.3 ±0.05 


21.8 ±0.06 


VVDS-2331 


24.0 


± 


0.09 


23.6 


± 


0.07 


23.4 


± 


0.07 


22.7 


± 


0.06 


24.0 


± 


0.03 


23.6 


± 


0.02 


23.2 


± 


0.02 


22.9 


± 


0.01 


22.3 


± 


0.03 


22.3 ±0.10* 


21.6 ±0.11 


VVDS-6913 




















21.7 


± 


0.14 


23.1 


± 


0.02 


22.6 


± 


0.01 


22.3 


± 


0.01 


22.1 


± 


0.01 


21.3 


± 


0.00 






VVDS-5726 




















22.4 




0.17 


24.3 


± 


0.04 


23.6 




0.02 


23.2 


± 


0.01 


22.7 


± 


0.01 


22.1 


± 


0.01 






VVDS-4252 


22.5 


± 


0.03 


22.4 


± 


0.02 


22.1 


± 


0.01 


22.0 


± 


0.17 


22.9 


± 


0.01 


22.4 


± 


0.01 


22.1 




0.01 


22.0 


± 


0.00 


21.5 


± 


0.00 


21.5 ±0.08* 


21.0 ±0.08* 


VVDS-4103 




















22.3 


± 


0.17 


23.2 


± 


0.02 


22.7 


± 


0.01 


22.4 


± 


0.01 


22.3 


± 


0.01 


22.1 


± 


0.01 






VVDS-4167 




















21.9 


± 


0.10 


































21.3 ±0.20* 


VVDS-7106 


22.9 


± 


0.04 


22.8 


± 


0.04 


22.7 




0.04 


22.5 


± 


0.04 


23.0 


± 


0.01 


22.7 


± 


0.01 


22.6 


± 


0.01 


22.5 




0.01 


22.4 


± 


0.02 


22.0 ±0.15* 


22.1 ±0.15* 


VVDS-6027 


23.6 




0.07 


23.2 


± 


0.05 


23.1 




0.05 


22.9 


± 


0.08 


24.0 


± 


0.04 


23.4 




0.02 


23.3 




0.02 


23.1 


± 


0.02 


22.8 


± 


0.05 


22.6 ±0.50* 


23.2 ±0.50* 


VVDS-1328 


24.3 


± 


0.07 


24.1 


± 


0.06 


24.2 




0.08 


23.9 


± 


0.12 


24.2 


± 


0.04 


23.8 


± 


0.03 


23.6 




0.03 


23.5 


± 


0.03 


23.6 


± 


0.10 


23.1 ±0.20 


22.2 ± 8.03 



magnitudes from the UKIDSS survey. 
Table 4. Stehar population properties 



Galaxy (1) L1500 (2) SFR^^^ (3) SFR,,d (4) E(B-V),ed (5) age (6) 



VVDS-1235 


6±2 


2.3 ±0.2 


1.8 ±0.2 


11 ±5 


0.28 ±0.12 


0.68 ±0.3 


VVDS-2331 


11 ±3 


3.7 ±0.1 


3 ±0.1 


27± 14 


0.32 ±0.11 


0.54 ±0.3 


VVDS-6913 


43± 16 


8.7 ±0.4 


7.1 ±0.4 


49 ±22 


0.28 ±0.10 


0.96 ±0.5 


VVDS-5726 


50± 17 


3.1 ±0.4 


2.4 ±0.4 


42 ±28 


0.37 ±0.13 


1.26 ±0.7 


VVDS-4252 


21 ±7 


11.2 ±0.3 


9.4 ±0.2 


73 ±26 


0.28 ±0.08 


0.40 ±0.2 


VVDS-4103 


16 ±4 


9.4 ±0.5 


7.7 ±0.4 


91 ±49 


0.28 ±0.03 


1.75 ± 1.6 


VVDS-4167 


25 ±8 


12.5 ±0.5 


10.6 ±0.4 


129 ± 46 


0.29 ±0.07 


0.29 ±0.1 


VVDS-7106 


8±2 


12.9 ±0.1 


10.6 ±0.1 


44± 11 


0.16 ±0.04 


0.27 ±0.1 


VVDS-6027 


11 ±4 


5.1 ±0.2 


4.1 ±0.2 


33 ±12 


0.24 ±0.08 


0.46 ±0.3 


VVDS-1328 


11 ±6 


4.4 ±0.2 


3.5 ±0.2 


14 ±6 


0.19 ±0.11 


0.88 ±0.6 



The columns are as follows: (1) Stellar mass (10^ Mq), (2) Luminosity at 1500A (xlO+2^ ergs s~^ Hz~i), (3) raw SFR derived from 
UV luminosity (Mq yr"-*^), (4) deredenned SFR derived from the full SED (averaged in the last 10^ years), (5) reddening suffered by 
the stars, (6) age of the oldest stehar populations (Gyr). 



filters from the CFHTLS T003 release ( |Ilbert et al]|2006 >. 
Finally, we have also used magnitudes in J and K bands 
from the UKIDSS survey (L awrence et al.||2007] . All the 
photometric data are presented in Table |3] 

Table [4] list the properties of the stellar population for 
the sample of the 10 galaxies derived from the SEDs. The 
stellar masses, the reddening suffered by stars, and the ages 
of the oldest stellar populations are estimated thanks to the 
full optical and near-infrared photometry available for our 
targets. The photometric points are compared to those of a 
library of synthetic stellar population spectra based on the 
latest 2007 models of Chariot k Bruzual (Bruzual 2007, lat- 



est models in preparation) models, using the Chabrier (2003 ) 
IMF. The synthetic stellar populations are based on expo- 
nential declining star formation histories with the addition 
of random secondary bursts. We calculated the probability 
distribution function as e*^~^ for each observed galaxy com- 
pared to all models in the library. We take as an estimate of 
the stellar mass the median of this distribution. The results 
are self-consistently corrected for the effects of age, metal- 
licity and dust. Degeneracies between these parameters are 
reflected in the error bars of the derived quantities. Fig. [2] 
shows the comparison between the observed SEDs and our 
best-fitting stellar population models. 



3.2.2 Star formation rates from photometry 

Ultraviolet-derived star formation rates were calculated 
from the i^-band photometry (see Table [sj. Given the red- 
shift of our targets, the i^-band is sensitive to the rest-frame 
UV continuum flux around 1500 A; in the absence of dust. 



the UV continuum from a young stellar population is ap- 
proximately flat if the flux is expressed in frequency units 
(fu)- The mean luminosity in frequency units (ergs s~^ Hz~^ 
) is calculated using the following equation: 



-.(-OA-m(AB)) 



D 



L(cm) 

l + z ■ 



(2) 



We can then deduce SFRuv using the following : 
SFR uv{Mq yr"^) = 0.83x10"^^ L1500 (ergs s"^ Hz-^).(3) 



This equation is the Kennicutt (1998) calibration renormal- 
ized from |Salpete7| ( |1955| ) to |Chabrier| ( [2QQ3t IMF. Table [i] 
list the value of the star formation rate computed from the 
UV rest-frame continuum emission. 

Additionally the star formation rates, averaged on the 
last hundred million years, SFRsed, can be also derived 
from the full SEDs using the same approach as for the stel- 
lar masses, reddening, and ages. They are also indicated in 
Table |4] for the 10 galaxies of our sample. 



3.3 Reddening and dust corrected star formation 
rate 

The UV star formation rates may be strongly affected by 
extinction, its comparison with the Ha star formation rate 
may thus provide an estimate of the amount of dust. As- 
suming that all of the ionizing photons are reprocessed into 
nebular lines, therefore SFR^^"^ = SFR'f^y after correc- 
tion for dust extinction, we thus computed the E{B — V) by 
solving the following system of two equations: 
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Figure 2. Rest-frame SEDs from the broad-band CFHT optical and near-infrared NTT/SOFI (or UKIDSS) photometry for our sample 
of 10 galaxies (see also Table [sj. The best-fitting stellar population model is overlaid. 
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SFR corr 
SFR corr 



SFRhc • IQ^-^^^^-^'^gas fc^(0.6563) 

SFRuv • lO^-'^^'^^-^has fc^(o.i5oo) 



(4) 



where the obscuration curve for the ionized gas, k^(A), is 
given by |Calzetti| (|2QQ1|) : 



r (A) = 2.659 (-2.156 + 1.509/A - 0.198/A^ 
+0.011/A2) +4.05 

0.12 fim<X< 0.63 fim, 
F(A) = 2.659 (-1.857 + 1.040/A) + 4.05 

0.63 fim<X< 2.20 fim. 



(5) 



We, then, derived the deredenned star formation rate, 
SFRcorr, hsted in Table [2] 



3.4 Gas mass and gas mass fraction 

The relation between the star formation surface density 
^SFR (in Mq yr~^ kpc~^) and the gas sur face density T ig as 



(in Mq pc ) is given by the Schmidt law (Kennicutt|l998 ) 



(6) 



which has been renormalized to the Chabrier (2003) IMF. 



We thus calculate the star formation surface density by di- 
viding the deredenned star formation rate (SFRcorr) by the 
projected area of the galaxy (see Table [s]). We used the two- 
dimensional Ha emission distribution, deconvolved with the 
point spread function (PSF), to directly estimate the spatial 
extend of the ionized gas, Agas, so that Mgas — Agas X YjgaS- 
We then calculate the gas surface density by invert- 
ing the above equation, and finally calculate the gas mass 
by multiplying the gas surface density with the projected 
area. We also derived the gas mass fraction, which is /i = 
Mgas /{Mgas + M^), where M* is the stellar masses deduced 
from the SED modelling. The gas mass and the gass mass 
fraction are given in Table [2] We found that six of our galax- 
ies have high values of the gas fraction (/i > 0.5). 



3.5 Virial mass 

While the global velocity dispersion aiD alone does not al- 
low us to distinguish between any systematic velocity shift 
across the galaxy and any random motions, it can provide 
a reasonable estimate of the mass of the object. Assuming 
a virialized motion around the morphological galaxy centre, 
we can use this am as a crude estimate of the mass within 
the largest radius rgas (estimated from the 2a width of the 
Ha flux profile along the major axis and deconvolved with 
the seeing, see Table [s]) for which a velocity is measured 
using the formula: 



Mvir{criD) - 



Ca^Y^Tgas 

G 



(7) 



(where C = 3.4 for of a gas-rich disk with an average incli- 
nation angle; see |Erb et al.| ( [2QQ6bt ). In Tablej?] we indicate 
dynamical masses inferred from the a id- The values span 
a range from ^-^1-26x10^° M©, with four galaxies having a 
virial mass > 1x10^^ 



4 PROPERTIES FROM SPATIALLY 
RESOLVED KINEMATICS 

4.1 Ha kinematics 

We produced two-dimensional maps of the dynamics for the 
galaxies in our sample by using E3D, the Euro3D visual- 
ization tool (Sanchez et al. 2004), and a code for the fit- 



ting and analysis of kinematics (e.g. Sanchez et al. 2004 



2005). Our main goal was to determine the kinematics of 
the ionised gas using the strongest emission line. Ha. For 
each spatial pixel (spaxel) we fit the Ha emission line pro- 
file with a single Gaussian function, in order to characterize 
the emission line, and a pedestal to characterise any spectral 
continuum. We first smoothed the reduced cubes spatially 
with a Gaussian of FWHM = 3 pixels (0'.'37). The line flux, 
FWHM, central wavelength and the continuum (pedestal) 
were then fitted. From the results of this fitting we obtained 
maps (see Fig. [s] and Fig. [3]) of the Ha emission line inten- 
sity, the relative radial velocity map {Vr) and the velocity 
dispersion (a). The dispersion map was corrected for the 
contribution of the instrumental dispersion, as determined 
from the FWHM of unblended and unresolved sky lines. Er- 
ror maps for the velocity and dispersion measurements were 
also computed. The errors are dominated by the effect of 
random noise in fitting the line profile which produces large 
errors at low S/N. Errors ranged from 4 to 25 km s~^ for 
both the radial velocity and velocity dispersion maps. 



4.2 Gas morphology 

Fig. [3] and Fig. [4] contain images of the Ha emission line in- 
tensity for each of the galaxies. We also plot the morphology 
of the stellar component mapping by the CFTHLS /-band 
continuum. Nine of the objects have a distinct single peak 
in the Ha distribution. This peak emission of the star for- 
mation from Ha is coincident (or near-coincident) with the 
peak stellar continuum emission in the / band. Only one 
of the targets, VVDS-1235, has distinct resolved knot-like 
structures in the Ha distribution which are not seen in the 
stellar emission. 

Given the similarity between the distribution of the Ha 
flux in our IFU observations and the continuum in the opti- 
cal images, we decided to use our line flux maps, which have 
the better resolution, to infer the morphological parameters: 
the centre (as defined by the position where the flux is max- 
imum); the position angle (PA), and the inclination. The 
inclination is calculated as cos~^(6/a), where a and h are 
the radius of major and minor axis, respectively, which are 
measured as the la profile of the flux map. The distance 
between the kinematical and the morphological centres is 
given in Table [5] We also measure a kinematical PA in case 
it is slightly offset from the morphological PA. The kine- 
matical PA can be measured in two consistent ways: or we 
look for the major axis as the direction where the velocity 
gradient is maximum, or we look for the minor axis as the 
direction where the velocity gradient is minimum. The dif- 
ference between the kinematical and morphological PAs is 
given m Table |5l 



8 M. Lemoine-Busserolle and F. Lamareille 

Table 5. Morphological parameters and initial guesses for the rotation modehing. 



Galaxy 


Ac (a) 


PA (b) 


APA^as (C) 


i(d) 


Tgas (e) 


Agas (f) 


VVDS-1235 




42°* 




70° 


9.8 ± 1.4 


63 


VVDS-2331 


2.4 


-69° 


-41° 


52° 


9.2± 1.6t 


120 


VVDS-6913 


8.5 


-99° 


-4° 


63° 


20.2 ±0.8 


348 


VVDS-5726 


1.7 


-18° 


-10° 


24° 


6.7 ± 1.0 


112 


VVDS-4252 


1.8 


45° 


-3° 


50° 


10.6 ± 1.1 


225 


VVDS-4103 


2.6 


-167° 


-18° 


53° 


13.4 ± 1.2 


315 


VVDS-4167 


2.8 


162° 


-10° 


62° 


16.6 ± 1.4 


468 


VVDS-7106 


0.0 


-119° 


0° 


37° 


7.7 ± 1.6 


142 


VVDS-6027 




175°* 




41° 


6.9 ± 1.1 


108 


VVDS-1328 


1.2 


-27° 


30° 


20° 


7.1 ± 1.1 


127 



The columns are as follows: (a) distance between the maximum of Ho; distribution and the kinematic centre (kpc); (b) initial guess 
kinematic position angle; (c) difference between the morphological (inferred from the Ho; distribution) and kinematic position angles; 
(d) inclination inferred from the Ho; distribution (0° for face-on); (e) 2sigma profile of ¥La distribution map along the major axis (kpc). 
Uncertainty represents half of the PSF correction; (f) projected area of the nebular emission (kpc^). 
*: for objects which do not show a rotation, we give the morphological PA. 

for VVDS-2331, we give the radius calculated with respect to the morphological major axis, not the kinematical major axis. 

Table 6. Best-fit parameters from the rotation modelling. 



Galaxy 




imod (b) 


cro (c) 


Vs (d) 


Vrot (e) 


rc 


(f) 


res(v) (g) 


res(d) (h) 


VVDS-2331 


-64° ± 21 


51° ±33 


38 ±41 


-8 ±21 


193 ± 72 


4.8 


±8 


-1±43 


11 ±34 


VVDS-6913 


-113° ±28 


62° ± 19 


44 ±22 


-6 ±19 


139 ± 30 


10.7 ±3 


-2 ±30 


20 ±23 


VVDS-5726 


-6° ± 22 


24° ±7 


60 ±35 


-17± 15 


292 ± 71 


1.5 


±2 


-2 ±25 


5± 19 


VVDS-4252 


32° ± 11 


48° ± 15 


100 ± 12 


-2 ±7 


130 ± 23 


4.5 


±2 


4± 10 


0± 11 


VVDS-4103 


167° ± 23 


79° ±5 


58 ±37 


-10± 12 


118 ±27 


8.5 


±3 


-2 ±20 


1±28 


VVDS-4167 


-176° ±16 


72° ± 23 


54 ±25 


-11 ± 12 


257 ± 34 


10.5 


±2 


2± 14 


5± 18 


VVDS-7106 


-134° ±35 


37° ± 26 


82 ±6 


7±6 


105 ± 71 


9.4 


±4 


2±7 


-2 ±10 


VVDS-1328 


-10° ±25 


20° ±7 


36 ±17 


17± 10 


182 ± 56 


0.5 


±4 


-1±11 


7± 17 



The columns are as follows: (a) position angle (b) inclination (c) local velocity dispersion (km s"-*^) (d) systemic velocity (km s"-*^) (e) 
maximum velocity (km s"-*^) (f) radius where the maximum velocity is reached (kpc) (g) mean value and rms of the residual velocity 
map (km s"-*^) (h) mean value and rms of the residual dispersion map (km s"-*^). 

In reality, the spatial resolution is limited by the see- 
ing and the spaxel size. The observed radial velocity is thus 
the weighted convolution of the true radial velocity by the 
PSF. This PSF is modelled as a 2D Gaussian and takes into 
account the 3-pixel spatial smoothing. The weights come 
from the flux map of the line used to compute the velocity 
map: a spaxel where the observed line flux is negligible will 
not contribute to the convolution. Additionally, the observed 
velocity dispersion accounts not only for the true physical 
dispersion, but also for the variations in the velocity field in- 
side the width of the PSF. We generate a map of the velocity 
dispersion where the velocity gradient (determined from the 
modelling) has been subtracted. This yields a better mea- 
sure of the intrinsic local velocity dispersion, rather than a 
raw value of a which may be inflated if the velocity gradi- 
ent across a spaxel is large. From this map, we compute the 
flux- weighted mean velocity dispersion, ao (Table [gJ. 

The beam smearing introduced by the PSF causes two 
significant effects. First, the observed velocity map appears 
smoothed, therefore the velocity gradient and the maximum 
velocity are underestimated. Secondly, the observed disper- 
sion map shows a peak of dispersion near to the kinematic 
centre. We can reproduce modelled velocity and dispersion 
maps by applying mathematically the same weighted convo- 
lution to the ideal velocity field. A better modelling of the 
beam smearing would be to convolve our modelled veloc- 



4.3 Rotation modelling 

Each galaxy is modelled by a pure, infinitely thin, rotating 
disk. The parameters of the model are the kinematic cen- 
tre (xo,2/o), the position angle (PA), the inclination (i), the 
velocity offset (Vs) of the centre relative to the integrated 
spectrum, and the velocity curve Vc{r) where r is the ra- 
dius from the kinematic centre. We also use the true phys- 
ical velocity dispersion ctq as a model parameter, which we 
assume is constant, and represents the thickness of the ro- 
tating disk. Below its maximum, the velocity is assumed to 
increase linearly with the radius and then remain constant 
( Wright et aL]|2QQ7| ). The radial velocity V for any pomt is 



then defined with standard projection equations. Note that 
the PA gives the direction of positive radial velocities. The 
velocity offset accounts for redshift uncertainties. We define 
the velocity curve by two parameters, following [Wright et al.| 



(2007): the maximum velocity Vmax and the radius Vc where 
this maximum velocity is achieved. In this model, below rc 
the velocity is defined as: 

Vc{r) = - X V^ax (8) 
rc 

The model computes Vrot, the asymptotic maximum rota- 
tion velocity at the plateau of the rotation curve, corrected 
for inclination (Vrot = Vrot/ sini). 
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Figure 3. From left to right: (a) /-band CFHT image, (b) Ho; flux map, (c) Ha velocity fleld, and (d) Ho; velocity dispersion map 
obtained from Gaussian flts to the SINFONI data cubes after smoothing spatially with a two-dimensional Gaussian of FWHM = 3 pixels. 
The /-band image and Ho; maps are color-coded with a linear scahng such that the values increase from light to dark. These data have 
been acquired with the 125x250mas sampling conflguration of SINFONI, in seeing limited mode. An angular size of l^^corresponds to 
~ 8 kpc at the redshifts of most of the objects. The direction of North-East, the VVDS number identiflcation and VIMOS-based redshift 
are indicated for each galaxy. 



ity and dispersion maps with the PSF in each wavelength 
channel, weighted by the intrinsic flux distributions of the 



systems at this wavelength (Emsellem, Monnet, & Bacon 
1994| . However, for such high redshift galaxies for which 



the intrinsic flux distributions is unknown and taking into 
account the quality of the data, it is not possible to model 
the velocity and dispersion maps to such a high precision. 
The modelled maps are compared to the observed maps by 
a x^— minimization. The parameters which minimize the 
are computed in two successive grids with increasing resolu- 



tion. First guesses for the kinematic centre, the PA, and the 
inclination are set from the morphology. 



4.4 Dynamical classes, masses and disk stability 

In order to characterize the galaxies in our sample, we used 
a kinematical classification based on the comparison of the 
morphology of the stars (showing by the CFHT i-band im- 
age) with the morphology of the Ha flux distribution. We 
classified the galaxies through analysis of the observed ve- 
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Figure 4. Similar as Fig[3]for five other galaxies. 



locity map, the observed velocity dispersion map (a, cor- 
rected from the instrumental effect), and the velocity resid- 
uals map. 

On the 10 galaxies of our sample, only two objects 
present a dynamical structure which is not compatible with 
rotation and which has not been successfully fitted by the 
simple rotating disk model. One object, VVDS-1235, has a 
complex kinematics with multiple velocity shear and sev- 
eral peak regions in the cr-map, and we defined as a merger 
and the other, VVDS-6027, shows negligible velocity shear 
(defined as featureless^ with the possibility of being a face- 
on disk). The remaining eight objects appear compatible 
with rotating disks and we used the correlation between the 
integrated line width, ctid, and the maximum velocity dis- 
persion estimated by the best-fitting rotating model to dis- 



tinguish rotation-dominated (RD) galaxies from dispersion- 
dominated (DD) galaxies ( Weiner et al.|200"6 ). Fig. [s] shows 
that four galaxies of our sample (VVDS-2331, VVDS-5726, 
VVDS-4167, and VVDS-1328) are RD disk (open circle) 
and four others (VVDS-4252, VVDS-4103, VVDS-7106 and 
VVDS-6913) are DD disks (filled symbols). Among the DD 
galaxies we highlighted using a filled square symbol, the 
galaxy VVDS-6913 which is consistent with being the relic 
of a major merging event. And using a filled triangle, the 
galaxy VVDS-7106 for which the error on the Vrot is ex- 
tremely high. 

We also infer the ratios of the maximum velocity over 
the local velocity dispersion Vrot/cfQ in order to characterize 
the dominance of the rotation versus the disordered motions 
of the gas. We determine the ratio rc/rgas between the ra- 
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Figure 6. Kinematic best-fit models of galaxies of our sample (left panel) for both the velocity and the dispersion, after adding the effect 
of the beam smearing. The difference for both the velocity and the dispersion between the observed map and the best-fit model map is 
also shown and may display non-negligible random motions in the gas. An one-dimensional rotation curve extracted using an 'idealized' 
slit from the observed Ho; velocity map, overlaid with the best-fitting model, is shown in the right panel. 



dius of the plateau, as estimated by the rotation modelling, 
and the radius of nebular emission detected in the SINFONI 
data to emphasize how close or far to the centre of the object 
the maximum of the velocity is reached. We diskuss this two 
ratios for each galaxy in section [s] The values for Vrot / cfq , 
Tc/rgas and the classification are shown in Table [t] 

For objects where the velocity shear is well fitted by the 
simple rotation model, we compute the dynamical masses 
assuming that the object is a circularly rotating disk. Thus 
the total mass within Vc for which the maximum rotational 
velocity Vrot is reached, is approximately described by: 

Mdyn = VLrc/G (9) 

where Vrot has been correction for inclination effects. Both 
this asymptotic velocity, and the radius of the turnover (rc) 
are inferred from our model fits to the observed velocity 
maps to correct for the the effect of beam-smearing (due to 
the seeing). The results are given in Table [7] 

We also computed the Toomre parameter 'Q' (Toomre 
1964) which is used to quantify the stability of the disk 
against gravitational collapse. Typically, values of Q > 1 
represent dynamically stable systems for which the total 
baryonic mass can be supported by the observed rotational 
motion alone, while values of Q < 1 suggest that the ob- 
served baryonic total mass exceeds that which can be sup- 
ported by the observed rotational velocity of the disk, in- 
dicating either that the disk itself is unstable or that non- 
rotational motions contribute a significant degree of support 



to the system. Finally, Q 1 represents the case of marginal 
stability. In order to estimate 'Q', we considered the ratio 
of masses within the radius probed by the SINFONI detec- 
tions (vdisk — rgas) and therefore computed the dynamical 
masses within the gas radius. We do not take the values 
of the dynamical masses computed from the rotationally- 
supported thin disk model (Mdyn)- The total baryonic mass 
is Mdisk = Mgas + M^. We thus used the following relation: 

Q = *™ = (10) 

CMdisk/rdisk ^(Mstar + Mgas)/rgas 

Results are given in Table [T] 



5 DYNAMICAL PROPERTIES AND NATURE 
OF INDIVIDUAL GALAXIES 

In this Section, we summarize the properties of each galaxy, 
derived from its spectral energy distribution, mass, star for- 
mation rate and dynamical structure. We use this to provide 
an insight into the nature and evolutionary status of each 
the 10 galaxies in our sample. 

5.1 VVDS-1235 

This galaxy has the lowest redshift of our sample and the 
lowest stellar mass 0.6 x 10^° Mq) . It exhibits complex 
kinematics, and it is impossible to find a clear centre for this 
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Figure 7. Similar as Fig[6]for four other galaxies. 



clumpy galaxy. Its cr-map shows several peaks which seem 
to match with the knots in the Ha intensity map. The peak 
with the maximum Ha flux shows a very high dispersion of 
240 km s~^. Our simple rotating disk model did not succeed 
in reproducing its complex velocity field; we interpret this 
galaxy as being a merging system. 



5.2 VVDS-2331 

This galaxy shows a velocity gradient of ±50 km s~^ in 
the south-west /north-east direction (position angle: —69°). 
However, we only cover the central part since the plateau is 
not reached in our SINFONI observation. Its cr-map show a 
peak at 120 km s~^ slightly shifted from the Ha flux peak. 
However this peak also corresponds to a region with high 
noise. The rotational disk model fits the observed inner gra- 
dient relatively well, and we found a plateau velocity, cor- 
rected from beam smearing, of Vrot 197 km s~^ within 
4.8 kpc. We found a virial mass of ^ 6 x 10^° Mq and a dy- 
namical mass of ^ 0.4 x 10^° M©, estimated from the Vrot 
inferred by the rotation modelling. It is one of the objects for 
which we found the lowest value of local velocity dispersion 
(< 40 km s"^), leading to a Vrot/cro ^ 5. VVDS-2331 also 
seems to have a large amount of gas (/i ^ 70%), a high SFR 
and a young stellar population. It is therefore a good candi- 
date for being a cold rotating disk (Q ~ 3), where rapid star 
formation is occurring due to the injection of large quantities 
of cold gas. 



5.3 VVDS-6913 



The Ha-map shows two components, one of which is coinci- 
dent with the peak of the Ha flux. However, the kinematical 
centre is located exactly between the two components, at ^ 
8 kpc from the centre of the main blob. It has a smooth ve- 
locity gradient and a large peak in the centre of its cr-map. 
The whole velocity map, including the two blobs, is well fit- 
ted by the rotation model (plateau velocity of ^ 139 km s~^ 
reached at 10.7 kpc and Vrot/cro ^ 3). However the small 
arms linked to the large peak in the cr-map reveal a per- 
turbed rotation. The peak of the dispersion and the arms 
are still visible in the residual map of the dispersion (see 
Fig. [g]), which emphasize additional random motions (up to 
70 km s~^ with a mean of 20 km s~^). This morphology 
can also be explained by the presence of a relic of a previous 
merger event, or by being a clumpy disk galaxy. VVDS-6913 
has the second highest virial mass 20 x 10^° ^o) and dy- 
namical mass ('-^ 5 X 10^° Mq) of the galaxies in our sample. 
It has a reasonable SFR and also a sizeable stellar mass 
4 X 10^° ^o) and the remaining quantity of gas (/x < 35%) 
indicates that it has already turned a large amount of its 
gas into stars. We classify it as a DD rotating disk, since its 
disk stability appears marginal (Q ^ 2). 
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Figure 5. Maximum velocity from rotation modeling (log Vrot) 
vs integrated line width (log ctid). open circle: RD rotating disks 
and filled symbols: DD rotating disks. Among the DD rotating 
disks we highlighted in filled triangle: VVDS-7106, and in filled 
square: VVDS-6913. The Hnes are the 1:1 line and the Rix et al. 
(1997) (7 = O-GVc line. 



5.4 VVDS-5726 

This nearly edge-on galaxy has its Ha-map and cr-map 
peaked at the centre. The general shape of the velocity gra- 
dient is well reproduced by the rotational disk model with a 
maximum velocity of Vrot ^ 292 km s~^ reached at 1.5 kpc 
and a local velocity dispersion of ^ 60 km s~^ {Vrot/cfQ 5). 
We found a virial mass of 9 x 10^° Mq and a dynamical 
mass of ^ 3 X 10^° Mq for this object. This galaxy has the 
highest stellar mass (^ 4 x 10^° M©), the lowest gas fraction 
(< 30%) and one of the oldest stellar population (1.26 Gyr) 
of all the objects in our sample. Also taking into account its 
moderate SFR, we believe that this galaxy has a rotating 
disk of cold gas (Q ^ 3) and it is at a later evolutionary 
state than most of the objects in our sample. 



5.5 VVDS-4252 

The resolved velocity structure present a smoothly varying 
shear along the major axis, including the beginning of a 
plateau at a radial velocity of ±54 km s~^. There is a very 
good agreement between the best-fit rotating disk model, 
with a maximum Vrot ^ 130 km s~^ reached at 4.5 kpc, and 
the observed rotation curve (see Fig.]?]). The residual maps 
show very low mean and rms values (see Table [6]). However, 
the velocity dispersion map is broadly peaked at the centre 
of the galaxy and the high value local velocity dispersion 
{r^ 100 km s~^) probably indicates that the galaxy does not 
have a dynamically cold rotating disk of ionized gas, due 
to significant disordered motion in the gas {Vrot/cfQ ^ 1). 
We also found that the virial mass {^^ 14 x 10^° M©) is 
much higher than the dynamical mass (^ 2 x 10^° M©). 



It has a young stellar population, a very high gas fraction 
(^ 87%) and seems to be undergoing a strong episode of 
star formation. VVDS-4252 is probably consistent with a 
rotation in a heated disk (DD unstable rotating disk (Q ~ 
0.4)). 

5.6 VVDS-4103 

This galaxy has a velocity gradient which can be fitted by a 
simple rotating disk model, and presents a clear peak in its 
cr-map which corresponds to the maximum of the Ha flux 
distribution. The best-fit disk model indicates a maximum 
velocity of Vrot 118 km s~^ km/s, reached at 8.5 kpc and 
a local velocity dispersion of ^ 58 km s~^. Hovewer, the 
peak in the residual dispersion map (see Fig. [t]) and the low 
ratio Vrot/cfQ ^ 2 indicates the presence of non-negligeable 
random motions in the gas. We also found that the virial 
mass {r^ 12 x 10^° Mq) is much higher than the dynamical 
mass ('-^ 3 X 10^° Mq). It has a extremely high SFR and gas 
fraction. VVDS-4103 must be experiencing a very strong 
burst of star formation, which might be the cause of non- 
negligible random motions of the gas peaked at the centre 
of the object. These galaxy properties are quite similar to 
the properties of VVDS-4252 and we also classify it as a DD 
unstable rotating disk (Q ^ 0.3). 

5.7 VVDS-4167 

This galaxy exhibits a smoothly varying shear along its ma- 
jor axis and has a clear peak in its cr-map corresponding to 
the maximum of the Ha flux distribution which is slightly 
shifted from the centre of the object. The velocity field is rel- 
atively well fitted by our simple rotational disk model with 
a maximum velocity of Vrot 257 km s~^ km/s, reached 
at 10.8 kpc, and a local velocity dispersion of 54 km s~^. 
This leads to a Vrot/cfQ 5 showing the dominance of the 
rotation over the random motions of the gas, even if the peak 
near the centre is also visible in the residual of the dispersion 
map (see Fig.]?]). We found a virial mass of 26 x 10^° Mq 
and a dynamical mass of - 17 X 10^° Mq for this object. 
It also has a sizeable and very young stellar population, a 
reasonable SFR and has turned slightly more than half of 
its gas into stars. We believe therefore that this galaxy can 
be consistent with a cold rotating disk (Q ^ 8). 

5.8 VVDS-7106 

This galaxy seems to have a low velocity shear (see Fig. [s]). 
As with VVDS-2331, we cover only the central part of this 
object and therefore the plateau has not been reached. The 
rotational disk model fits the observed inner gradient rel- 
atively well and gives an estimated maximum velocity of 
Vrot ^ 105 km s~^ km/s, reached at 9.4 kpc, and a local 
velocity dispersion of ^ 82 km s~^. We also found that the 
virial mass ('-^ 4 x 10^° Mq is similar to the dynamical mass 
(-- 2 X 10^° Mo), and that it has a ratio of > 1. We 
measure a low ratio of Vrot/cro ^ 1, showing the dominance 
of random motions in the gas as also indicted by the high 
local velocity dispersion. VVDS-7106 has the lowest stellar 
mass in our sample and an extremely young stellar popu- 
lation. More than 60% of its gas has been turned into star 
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and it has an average SFR. Taking into account all of its 
properties, we classify it as a DD rotating disk despite its 
marginal stable disk {Q ^ 1). 

5.9 VVDS-6027 

This object has a young stellar population with an average 
stellar mass of - 1 X 10^° Mo, a low gas fraction of (< 35%) 
and a small SFR. It presents an Ha flux distribution with 
a peak located at its centre and a low surface brightness 
emission tail seen also in the stellar population image (see 
Fig. [4]). This galaxy is relatively well-resolved but shows no 
strong evidence for spatially resolved velocity structure. We 
classify it as featureless (perhaps consistent with being a 
face-on galaxy). 



5.10 VVDS-1328 

As for VVDS-5726, this nearly edge-on galaxy has a ve- 
locity gradient which is well reproduced by the rotational 
disk model. However, the peak seen in the cr-map is slightly 
shifted from the peak in the Ha flux map (see Fig. [s]). 
The best-fit disk model indicates a maximum velocity of 
Vrot ^ 182 km s~^ km/s, reached at 0.5 kpc, and a local 
velocity dispersion of ^ 36 km s~^ (Vrot/cro ^ 5). We found 
a virial mass of 3 x 10^° Mq and a dynamical mass of ^ 
0.4 X 10^° Mq for this object. This galaxy has an average 
stellar mass of - 1 X 10^° Mo, a low gas fraction of (< 35%) 
and a small SFR. We believe that this 'mature' galaxy has 
a rotating disk of cold gas (Q 6). 



6 DISKUSSION 

Here, we study and diskuss the relations between the kine- 
matical and physical properties of the ionized gas and the 
physical properties of the stellar population for the 10 galax- 
ies of our sample. 

6.1 Relations between the recent-to-past star 
formation rate and other properties 

We calculated the recent-to-past star formation rate ratio 
{SFRcorr / SFRsed)^ which is the ratio between the star for- 
mation rate estimated from the ionized gas (using the Ha 
emission) and the one estimated from the best-model SED 
fitting (both of which have been dereddened). SFRsed is av- 
eraged on a timescale ten times longer than the SFRcorr, 
which is an indication of the instantaneous star forma- 
tion. Given these different timescales, the ratio between the 
two estimates gives information on the presence of recent 
burst of star formation above the standard continuous de- 
clining star formation. We investigated, thus, the correla- 
tion between the recent-to-past star formation rate ratio 
{SFRcorr/ SFRsed) and the gas mass fraction and the stellar 
mass. 



6.1.1 Relations between the recent-to-past star formation 
rate ratio and the gas mass fraction 

In Fig. [8] we plot the recent-to-past SFR ratio vs. the 
gas mass fraction for the 10 galaxies of our sample. We 
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Figure 8. Relation between the recent-to-past star formation 
rate ratio and the gas mass fraction. The symbols give the dy- 
namical classification as follows: open circle: RD rotating disks; 
filled symbols: DD rotating disks, with in particular filled trian- 
gle: VVDS-7106, and filled square: VVDS-6913; star: featureless. 
We found two classes (see text) and a least-square linear fit is 
performed for each of them. 



found that only galaxies with a significant amount of gas 
(/i '^>0.6), are undergoing episodes of strong star formation 
(bursts). From Fig. [s] a correlation between the recent-to- 
past SFR ratio and the gas mass fraction, given that we 
divide the sample into two classes is therefore clearly seen. 
The first one include the three objects which have the small- 
est recent-to-past SFR ratio. These three objects are: the 
face-on/featureless VVDS-6027 (star symbol), the DD ro- 
tating galaxy VVDS-7106 (filled triangle), and the RD ro- 
tating galaxy VVDS-4167 (open cirle), and therefore are not 
showing evidence of bursts of star formation. This first class 
seems to follow the relation (least square fit): 



logioiSFR corr /SFRs 



1.38(±0.45) X /i - 1.40(±0.24), 



and the second class, all the remaining objects, seems to 
follow the relation: 

logio{SFRcorr/ SFRsed) = 1.77(=b0.01) X /i - 1.03(±0.01). 

Can these two classes be related to underlying kinematical 
properties ? We find indeed that the ratio of the kinematical 
radius to the gas radius {rc/rgas, see Tabled has a mean 
value of 0.40 ± 0.01 for the galaxies with a high recent-to- 
past SFR ratio, and 0.94 ± 0.01 for the three others of the 
first class. It appears, thus, that the rather intuitive relation 
between the amount of recent star formation and the amount 
of ionized gas depends strongly on kinematics. 'Slow rotator' 
(intermediate redshift galaxies with a rc/rgas > 0.6) appear 
less efficient in converting a high mass fraction of ionized gas 
into a burst of star formation. This will have to be confirmed 
and refined with better data and more objects. 
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Table 7. Kinematical properties 



Galaxy 


shear 


Vrot 1 (^0 


M^ir (1) 


Mdyr. (1) 


f^gas 


Q(2) 


dyn. class (3) 


VVDS-1235 


92±13 




4±2 








merger. 


VVDS-2331 


99±17 


5.1 ±2.9 


6.4 ±2.9 


0.4 ±0.7 


0.52 ±0.87 


3.3 ±0.9 


RD 


VVDS-6913 


139±7 


3.2 ±0.8 


20 ±6.5 


5±3 


0.53 ±0.15 


2.4 ±0.6 


DD 


VVDS-5726 


95±9 


4.9 ±0.5 


8.7 ±0.3 


3±4 


0.22 ±0.3 


3.3 ±0.5 


RD 


VVDS-4252 


54±9 


1.3 ±0.1 


14 ±8.7 


2± 1 


0.42 ±0.19 


0.4 ±0.7 


DD 


VVDS-4103 


83±16 


2.0 ±0.7 


12± 13 


3±2 


0.63 ±0.23 


0.3 ±0.8 


DD 


VVDS-4167 


140±8 


4.8 ±0.8 


26 ±2.2 


17 ±5 


0.65 ±0.13 


7.5 ±0.5 


RD 


VVDS-7106 


12±14 


1.3 ±0.4 


4.1 ±0.7 


2±3 


1.22 ±0.58 


1.4 ± 1.0 


DD 


VVDS-6027 


<10±11 




1 ±0.5 








f.l. 


VVDS-1328 


58±8 


5.1 ±0.7 


3.4 ± 1.0 


0.4 ±3 


0.07 ±0.56 


5.7 ±0.7 


RD 



The columns are as follows: (1) 10^^ (2) Toomre parameter (Toomre 1964); (3) dynamical classification: RD for 
rotation- dominated rotating disk, DD for dispersion- dominated rotating disk, and f.l. for featureless. 




and all RD galaxies, except VVDS-2331, seems to follow the 
relation (least squares fit) 



Figure 9. Relation between the recent-to-past star formation 
rate ratio and the logarithm of the stellar mass. The symbols are 
the same than in Fig.js] 



6.1.2 Relations between the recent-to-past star formation 
rate ratio and the stellar mass 

In Fig. [9] we plot the recent-to-past SFR ratio vs. the stellar 
mass for the 10 galaxies of our sample. We found that the 
galaxies of our sample need higher stellar masses to reach 
higher recent star formation rates. The recent star formation 
to stellar mass efficiency is moreover significantly stronger 
for half of our sample. From Fig. [9] we see two possible corre- 
lations between the recent-to-past SFR ratio and the stellar 
mass. The first class includes objects dominated by recent 
star formation {SFRcorr / SFRsed > 1) - the merger object 
and all DD galaxies, except VVDS-6913 with the following 
relation (least squares fit): 

logioiSFRcorr/ SFRsed) = 0.072(±0.05) x M^IO^M©) 

-0.72(±0.4) 

The second class includes objects dominated by past star 
formation {SFRcorr/ SFRsed < 1): the featureless object 



logio{SFR 

corr /SFRs 



: 0.014(±0.03) X M^lO^Mo) 
-0.97(±0.7) 



The star formation to stellar mass efficiency is moreover 
stronger for objects experiencing a burst of recent star for- 
mation, and appears also to be related to the dominance of 
random motions in the gas. We indeed notice that the ratio 
of the maximum velocity to the velocity dispersion {Vrot /(Jo] 
see Table [7]) is different for the rotating disks galaxies inside 
the two classes: it has a mean value of 2.42 ± 0.5 for the first 
class, and 4.5±0.3 for the second one. We thus propose to re- 
fine our classification: disk galaxies of the first family would 
be called "hot /perturbed rotators", while the disk galaxies 
of the second family would be 'cold rotators'. Only hot ro- 
tators and objects dominate by random motions, like the 
merger seem to be able to experience a recent burst. The 
separation between cold rotators and hot/perturbed rota- 
tors would then be Vrot/o-Q ^ 4. This also will have to be 
confirmed and refined with better data and more objects. 



6.2 The stellar mass Tully-Fisher relation 

In Fig. [To] we show the Tully-Fisher relations for the 
1.2 < z < 1.5 galaxies showing evidence of rotation. In the top 
panel, we plotted the stellar mass versus the maximum ve- 
locity. As can be seen, there is a clear separation between the 
RD rotating disks (open circle) from the DD rotating disks 
(filled symbols). However, both groups of objects suggest a 
good correlation between the stellar mass and the rotation 
maximum speed similar to that seen in the local universe 
(sold line). By assuming the local slope of the Tully-Fisher 
relation from ( |Bell de Jong|2001 ), we found a best fitting 
zero point of -0.29±0.25 for the RD galaxies and 0.90±0.43 
for the DD galaxies. 

In the bottom panel of Fig. [lO] we replace the max- 
imum velocity by the S0.5 index {0.5Vrot + ^^0)°"^), which 
allows the disordered motions in the gas to be taken into 
account (Kassin et al. 2007). The separation between the 



RD rotating disks and the DD rotating disks is less obvious 
while adding the disordered motions. It seems that we also 
detect an evolution of the relation since Z'-^l (which is shown 
as the solid line in Fig. 10). 
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Figure 10. The stellar mass Tully-Fisher relation at 1.2<z<1.5, 
using the maximum velocity (top panel) and using the S'0.5 index 
(bottom panel; {Q.bV^^^ +'^n)^'^)- ^^P P^tnel, the solid line 



is the local relation from ( Bell de Jong||2001|>, while i n the 
bottom panel it is the z~l relation from f Kassin et al.|2007| >. The 
symbols are open cirle: RD rotating disks and filled symbols: DD 
rotating disks with in particular filled triangle: VVDS-7106, and 
filled square: VVDS-6913. 



6.3 Comparison with theoretical studies, with 
results from other intermediate- and 
high-redshift samples 

Previous studies at z > 2 have shown the presence of large 
quantity of gas in comparison to the stellar population in 
high-redshift galaxies ( Law et aL]|2009 Lemoine-Busserolle 
et al.|[2009). We found that six out of then galaxies of our 



Figure 11. Combined velocity (.So. 5 ) versus integrated hnewidth 
velocity dispersion. The diagonal lines are the same as in jWrighT] 
|et al.| ( |2007l >, i.e. a 1 : 1 line and the |Rix et ar] ( |1997t a = 0.6Vc 
line (Vc is the circular velocity). Linewidth auj and S'0.5 are 
correlated; the 0.5 pre-factor makes the combined velocity width 
a better estimate of velocity dispersion, so that the correlation 
is tighter and the galaxies closer to the 1 : 1 line. The sym- 
bols are open cirle: RD rotating disks and filled symbols: DD 
rotating disks with in particular filled triangle: VVDS-7106, and 
filled square: VVDS-6913; asterisk: 2; ~ 3 galaxies from |Lemoine-| 
[Busserolle et al.| ( |2009] > 



v/a for our sample is similar to the one found for the SINS 
galaxies ( 



et al. 



Forster Schreiber et al. (2006); Genzel 



et al. (2006 2008)), but more higher than the values found 
by Lemoine-Busserolle et al. (2009) (~ 0.4 — 1.5) and Law 



et al. 



2009) galaxies at z ^ 2 0.8). Lemoine-Busserolle 
39) found in their sample, 3 galaxies showing the 
presence of rotation. Fig[TT]plots the combined velocity scale 
5'o.5 against the ID linewidth am for our sample together 
with this three galaxies at 2; 3 — 4. It should give an idea 
of the evolution of the dynamical mass which take also into 
account the contribution from the random motions in the 
gas. Fig [11] shows an increase of the dynamical mass for 
z - 3 - 4 to z - 1.3 - 1.5 for the DD disks. 

Despite the low resolution of the morphological data 
(CFTHLS I-band images) which shows that, except VVDS- 
1235, all the galaxies seems to be disk like objects with a 
central bulge, the IFU data obtained on the sample of the 
10 star-forming galaxies at 1.0<z<1.5 have confirmed the 
presence of a velocity gradient resembling that expected for 
an internally rotating system, but has also reveals the pres- 
ence of non-negligible random motions in the gas. Recent 
studies at the same redshift have found similar results. [Bour-I 



sample have a gas fraction higher than 0.50. The typical 



naud et al. (2008) have study a clumpy galaxy at 2; = 1.6, 
which was previously classified as an on-going merger. Their 
study reveals that the kinematics properties (the large-scale 
rotation with Vrot 100 km s~^), the stellar properties 
(young age, stellar mass range) and the properties of the 
gas are very similar to those of the galaxies of our sample 
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(in particular VVDS-6913). Their numerical model shows 
that a such turbulent rotation disk, results more from the 
evolution of an unstable gas-rich disk galaxy than from a 



merging event. ( , Wright et al. 
^0081 



2007 



2009 



Bournaud et al. 



Five of the rotating disk of our sample have a Q > 2, 
which indicates the stability of their disk. Can this large 
stable disk systems evolve in present-day thin disks ? We 
already note that the disk galaxies in our sample (except 
the two galaxies VVDS-1328 and VVDS-5726 which already 
have similar properties to that of the present-day local spi- 
rals) possess a large disk of gas and younger stars (age < 1 
Gyr), which can further grow in mass by continued accre- 
tion. The lack of rotation observed in this objects in com- 
parison to the thin local disk is nowadays well known and 
taken into recent plausible formation and evolution of high- 
redshift thick disk models. These scenarios can predict that 
on the largest scales a velocity gradient tracing a rotation 
can be observed with high velocity dispersion. They also 
shown that before z ^ 2 the hot mode of accretion domi- 
nates, but after z ^ 2, the galaxy has a thin and extended 
disk component with Q 1.5 — 2, which indicates that the 
disk is marginally stable. They indicate also that after these 
earlier stages the galaxy enters a slow accretion phase and 
the disk evolves quiescently until z = ( Keres et al.||2005 



Dekel Birnboim|[2QQ6l |Ocvirk et aT]|2008| |Dekel et al.| 
2009a). The results presented in this paper support the hy- 



pothesis that stable gas-rich disks seen at intermediate and 
high-redshifts may internally evolve into present-day spirals. 



7 CONCLUSIONS 

We have presented the 2D kinematics and the physical 
properties of a sample of ten star-forming disk galaxies 
at 1.0<z<1.5. Among these objects, three (VVDS-4252, 
VVDS-4103 and VVDS-4167) are undergoing a strong burst 
of star formation. We found mainly four kinematical types 
in our sample which are: one merger (VVDS-1235); one fea- 
tureless (or face-on) galaxy (VVDS-6027); four DD rotating 
disks - VVDS-4252, VVDS-7106, VVDS-4103 and VVDS- 
6913 which is also consistent of being the relic on a ma- 
jor merging event (to which the disk system seems to have 
survived) or a good candidate for being a clumpy galaxy; 
and finally four other galaxies are RD disks — VVDS-2331, 
VVDS-4167, including VVDS-1328 and VVDS-5726 which 
are pure rotationally supported disks. These two rotating 
disks achieve a maximum velocity of ^ 180-290 km s~^ 
km/s within ~ 0.5-1 kpc, similar to local spirals with thin 
disk ( Sofue et al.|1999 ). Regarding most of the DD rotating 
disk, they display a plateau velocity range of 105-257 km 
s~^ km/s, certainly underestimated due to beam smearing. 
However, their plateau radii (4.5-10.8 kpc) derived from our 
rotating disk model are significantly higher than those de- 
rived for pure rotating disks and local spiral galaxies. We did 
not find any trace of AGN in the 10 objects. The galaxies of 
our sample have a relatively young stellar population (< 1 
Gyr) and possess a range of stellar mass of 0.6-5 10^° Mq. In 
addition, most of them have not yet converted the majority 
of their gas into stars (six galaxies have their gas fraction > 
50 per cent). Therefore, the galaxies which already have a 
stable disk (six of them) would have their final stellar mass 



similar to the present-day spirals, to which these rotating 
systems can be seen as precursors. We also investigated the 
stellar mass Tully-Fisher relation at 1.2<z<1.5 and found a 
scatter between the RD disk galaxies and the DD disk galax- 
ies. We also found for this two groups of objects a change 
of the zero point in comparison with the stellar mass Tully- 
Fisher relation in the local universe, but this is speculative 
considering the statistics of our sample. 

Although, we have presented results on the dynamical 
structure and the physical properties of a small sample of 
star-forming galaxies at 1.0<z<1.5, we have been able to 
investigate the dynamical type and the physical properties 
of late-type objects at this intermediate redshifts. Increas- 
ing the samples at different redshift ranges for which near- 
infrared IFS data and multi-wavelength broad-band pho- 
tometric data can be obtained, will definitely develop our 
understanding of the dynamical characteristics and nature 
of various systems (mergers, rotating disks, etc.) and there- 
fore will lead us to understand how galaxies have evolved to 
match the present-day Hubble sequence. 
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